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I. INTRODUCTION 



A. MID-ATLANTIC BIGHT OBJECTIVES 

Better knowledge of operating environments may allow any Navy to improve the 
use of its USW systems in different tactical scenarios. Even though undersea sound and its 
practical uses have been continuously explored since World War II, testing of other 
techniques are necessary in order to obtain important sea environment information which 
may help in the development of accurate Naval USW systems. To face this challenge, the 
Mid-Atlantic Bight (MAB) experiment was conducted. The MAB experiment was an 
intensive field study of the shelfbreak south of Long Island during the strongly stratified 
summer season whose basic objectives were: 

1. To provide detailed time series of the oceanographic conditions along typical 
major acoustic paths by using moored acoustical arrays. 

2. To provide a diversity of frequencies and source locations from broadband SUS 
sources. 

3. To determine the temperature structure of the ocean environment by using 
acoustic tomography. 

This study may provide the Navy an alternative and simple means of determining 
characteristics in unknown or poorly known operating environments. (Pickart et al., 1996) 

B. STUDY OBJECTIVES 

The objectives for reaching the fundamental goal of this study are as follows: 

1. To determine the bottom sound speed by measuring the coherence of two 
vertically displaced sensors in the shallow water channel and comparing this 
function with one obtained from a simple model equation over a range of 
frequencies. 

2. To determine the bottom density by comparing mode travel times of computed 
spectrograms for a phone from SUS explosive data and those obtained by 
employing the concept of an effective depth for an ideal waveguide. 



1 



3. To determine the bottom attenuation by finding mode attenuation as a function 
of range extracted from spectrograms and relating it to loss per bottom bounce. 
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II. EXPERIMENTAL/ENVIRONMENTAL DESCRIPTION 



The summer Mid- Atlantic Bight experiment consisted of a cruise from 19 July to 9 
August 1996 aboard the R/V Endeavor from the University of Rhode Island. The cruise 
included three major components: acoustic tomography, SEASOAR operations, and 
hydrography/tracers. An AXBT flight was also conducted during this time period, as well 
as a SUS drop fly-over. The acoustic portion of the experiment consisted of three 
components: 

1 . The moored acoustic transmission array, or tomography array 

2. The moored physical oceanography array, which provided detailed time series 
of the oceanographic conditions along our major acoustic path 

3. The broadband, explosive (MK61 SUS shot) source component. 

Multiple lines of SUS drops were deployed in the area by P-3 aircraft and monitored by two 
vertical line arrays (VLAs) moored on the shelf in ~85m water. (Pickart et al., 1996) 

To achieve the objectives of this thesis, eleven broadband SUS shot transmissions 
recorded at the NW VLA were selected. All of these were dropped along a single track 
from a point at latitude 40°16.2’N, longitude 71°10.2’W to a point at latitude 40°16.8’N, 
longitude 70°57.0’W. The sampling rate was 3906.25 Hz, the duration of one record was 
8.3886 sec and the number of samples/record was 32768. The arrangement and the location 
of the acoustical elements used for this study are shown in Figure 1 . This region of the 
experimental site has a nominal water depth of ~ 90 m with little variation (<10m) from the 
VLA to the positions of the selected SUS drops. The receiving VLA had 16 hydrophones 
of which only the upper 8 successfully recorded the SUS transmissions with an upper depth 
of 30.5 m and inter-element spacing of 3.5m. A typical sound speed profile measured near 
the trial site and used in this analysis is shown in Figure 2. 
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Figure 1. Arrangement and Location of Acoustic Elements Used at the Experimental Site. 
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Figure 2. Sound Speed Profile Measured Near the Trial Site and Used in This Study. 
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III. ANALYSIS METHODS 



A. INVERSION OF BOTTOM SOUND SPEED 

In determining the bottom sound speed a method was perfomied to determine the 
critical angle of the seabed from the vertical directionality of the ambient noise in the water 
column (Buckingham and Jones, 1987). The properties of the ambient noise that we are 
specifically concerned with are the vertical directionality and the spatial coherence of the 
fluctuations of two sensors vertically displaced in the noise field. These two quantities are 
related through a Fourier transform relationship, which is given by 

r (n)=!jv(0) exp[/Q cos(#)] sin(<9)(i 0 (1) 



where 0 is the polar angle measured down from the zenith and F(0) is a dimensionless 
directional density function which represents the distribution of the noise power in the 
vertical. F(0) can be normalized according to the condition 

|jV(0)sin(0)de = l. (2) 

This indicates that in the limit as Q — » 0 the value of the coherence function goes to unity. 

Buckingham and Jones provided a simple example for isotropic noise in which F(0) 
is independent of 0. According to Eqs. (1) and (2) the coherence function in this case is 
shown to be 



r isotr( Q ) = 



sin(Q) 

Q 



where Q is a dimensionless angular frequency defined as 

col 

Q = — , 



(3) 

(4) 



(0 is the angular frequency, c w is the sound speed in the water column, and / is the 
separation between the two sensor positions. This example is shown in Figure 3. 
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Figure 3. Example of Ideal Vertical Coherence Function Versus Frequency. 

This function predicts that the coherence of the noise at two sensors is real due to 
symmetry of the noise field about the horizontal. For an asymmetrical noise distribution in 
the vertical, the coherence function is complex with a large imaginary component. A 
symmetrical noise distribution is also predicted by Buckingham and Jones’ theoretical 
model of shallow water noise, where the coherence function for shallow water noise is an 
oscillatory function that decays with increasing . At a given frequency, the directional 
density function of the noise predicted by the model is closely approximated by 



where a x is independent of 9. Then, the noise propagates in directions close to the 

horizontal, with no acoustic energy traveling with grazing angles greater than 6 C (the critical 
grazing angle). The contribution of local, wind-driven sources to the overall noise field is 
approximately isotropic, that is, the directional density function of the continuous 
component of the noise can be expressed as 




( 5 ) 



0, otherwise 



F^(0) = a 2 , for 0 < 9 <7i , 



( 6 ) 
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where a 2 is independent of 9. By superposing Eqs. (5) and (6) we have a representative 
directional density function for the total noise field in shallow water: 



(a, +tf 2 )’ f° r 



3H ses (r 0N 



( 7 ) 



a ,, 



otherwise. 



Figure 4 shows a sketch of the vertical distribution expected in a shallow water 
environment. Through the normalization condition in Eq. (2), the critical grazing angle can 
be expressed in terms of the other two parameters as 

sin(e f ) = ^, (8) 

a 2 

and the model coherence function of the noise represented by the distribution in Eq. (7) is 
found from Eq. (1) to be 



/ \ a-, / \ 1 — a, 

rje) = -si„(a) + ^ 



sin^Qsin(0 c )j . 



Using this model, the spatial coherence of two vertically displaced sensors 
NW VLA was computed from extracted ambient noise data. By spatial coherence is 



( 9 ) 

of the 
meant 




Figure 4. Model of the Vertical Distribution of Ambient Noise in a Pekeris Waveguide. 
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the time-averaged crosscorrelation coefficient of the noise observed by hydrophones 
separated vertically in the water column, i.e., 



pM 

JpMpM ’ 



( 10 ) 



where P x> ,(co) represents the cross spectral density as a function of frequency for two 
separated phones and P xl (co) and P yy (®) represent the spectral density for each one of the 
separated phones. To search the two unknown parameters a 2 and 0 C in Eq. (9), an 
optimization procedure was implemented by finding a least squares fit between the real part 
of the measured coherence function, T meas , and the theoretical coherence function, T w . By 
defining the ambiguity function 

° = Z l r „ ( n >°2 • 6 . ) - Re al {r„ m .( n )}| 2 . (ii) 

n=l 

values for a 2 and Q c are obtained when O is a minimum. 

B. INVERSION OF BOTTOM DENSITY 

To determine the bottom density, the concept of an effective depth for an ideal 
waveguide was employed (Chapman and Ward, 1989). By considering the phase change of 
a plane wave reflected at a fluid/fluid boundary, Weston introduced the concept of the 
effective boundary depth, arguing that one could view the reflection as taking place at an 
imaginary pressure - release boundary located at a specific depth below the true boundary. 
He proposed that the normal modes of the Pekeris model should be similar to those of an 
ideal pressure - release waveguide whose depth was equal to the true water depth plus the 
effective boundary depth. The normal mode wavenumbers for the Pekeris model can then 
be estimated easily without numerical iteration. 

The Pekeris model consists of a uniform water layer of constant depth bounded 
above by a pressure - release surface and below by a uniform fluid half space having 
different density and sound speed. Consider a plane wave of frequency fin a fluid medium 
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of sound speed c w incident upon a plane boundary with another acoustic medium at angle 9 
relative to grazing, as depicted in Figure 5. 




The ray segments s and s ' are given by 



and 



s = 



AH 
sin 9 



2 A# cos 2 9 
sin 9 

and the phase of the reflected plane wave at points A and B is given by 
and 

<h = 2* w 5 + 4", 



( 12 ) 



( 13 ) 



( 14 ) 

( 15 ) 



in which k w 



(0 

— . Requiring <{>,=<()„ the effective depth of the boundary is defined by 
c w 






T , - i r 
2k w sin 9 



( 16 ) 
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For an isospeed water column with sound speed c w and density p overlying a 
homogeneous solid medium with compressional sound speed c ft , shear speed c s , and 
density p b , Chapman and Ward (1984), have shown (assuming c s < c M and bottom 
attenuation is negligible) that the effective depth is 

f-sinG f , ( 17 ) 

where H is the water depth and 0 f is the grazing critical angle. If c s « c w then Eq. (17) is 
reduced to 







fc } 


2 “ 


= H + 


1-2 


C s 













H=H+- 



g 



k,., sin6. 



( 18 ) 



The boundary conditions for this ideal waveguide require the pressure to vanish at 
the upper and lower interfaces. Correspondingly, eigenfunctions, \|/(z) , must satisfy 
(Medwin and Clay, 1998) 



YU)Uo= <t> and v|/(z)| z=// = 0 . 
Solutions of the eigenfunctions for these boundary conditions are 

vj/(z) = sin 



rmz 






( 19 ) 



( 20 ) 



and the vertical component of the wavenumber, y, must then satisfy the quantization 
condition 



mt 




( 21 ) 



where n is the integer and designates the mode number. The horizontal values of the 
wavenumber K n are then given by 




( 22 ) 
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and the requirement that K n be real gives the mode cut-off. By substituting y n from Eq. 



(21) and k w - — , K n can be written as 
c... 



K 1 



(- 



>2 ( \ 2 

nn 



H„ 



( 23 ) 



Furthermore, the group velocity of the n* mode is given by 

dc o 



= 



dK.. 



( 24 ) 



therefore. 



= • 



£i 1 


rl 


2 




0) ]j 


K C W J 




UeJ 



= c.. 



1- 



( V 

nnc,. 



( 25 ) 






The group velocity u n and the mode travel times are related by 

R 

K ~ > 

u 



( 26 ) 



where R is the horizontal distance between source and receiver and t n . is the travel time of 
the n* mode. Combining Eqs. (25) in (26) the mode travel times can be written 



t — 



1 


F 


( „„ \ 

nc w 





( 27 ) 



From the definition of effective depth, this becomes 



t.i/h-u 

c... 



«7tsin0. 



2nfH sin0 f 



+ g 



- 1/2 



( 28 ) 



which is the final expression for the mode travel times . Since we know the real depth H, 
and c w and 0 r were determined from the previous analysis, then we can invert for g by 
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measuring the mode travel times over a range of frequencies for different SUS 
transmissions. Assuming p w * 1, then g « p h . 

C. INVERSION OF BOTTOM ATTENUATION 

Mode attenuation as a function of range can be estimated from several computed 
SUS spectrograms. The mode attenuation per km, a„ , satisfies the relation 

ATL n =a n AR, ( 29 ) 



where ATL n and A R are the mode transmission loss due to attenuation and range between 
SUS drops, respectively. Neglecting volume attenuation in the water column, the decrease 
in mode amplitude is expected to be caused only by bottom attenuation and cylindrical 
spreading ( TL = lOlogR). Since mode n can be thought of as due to rays propagating at 
angle G n , the mode attenuation can be considered the result of multiple bottom reflection 
losses. 



The number of bottom bounces over 1 km of horizontal propagation can be found 
by considering the angle of elevation or depression 0 of the local direction of propagation 
of the waveform. Referring to Figure 6, the angle 0 and the diagonal distance for one 
bounce L , are related by 

2 LJ 

sin0 = , (30) 

A 

where H is the water depth. Defining the total diagonal distance L corresponding to R = 1 
km of horizontal propagation by 



L = 



R 

COS0 



( 31 ) 



then L and Z,, can be related in order to define the number of bottom bounces per km as 

N = — = R — Q . (32) 

L 2 H 



The attenuation for mode n can then be expected to be due to N(d n ) bottom bounces per 
km producing 

o. = w(e».). (M) 
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where 0„ is the corresponding propagation angle for mode n and a M (0 n )is the bottom loss 
per bounce for a single reflection at this angle. 




The reflection coefficient at a fluid-fluid interface for a plane wave traveling at 
angle 0 is defined by 

^p.c.sine + ipA, (34) 

Pb c b smQ- ip w b 2 

where c w and p w are the sound speed and density, respectively, in the upper medium, c b and 
P b are the sound speed and density, respectively, in the lower medium, and 



b 2 = 



\ 



cos 2 0 - 1 



vc u ,y 



1/2 



( 35 ) 
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If there is no attenuation, then |R| = 1 for incidence beyond the critical angle. However, if 
the bottom is lossy then 

|ft(ej=i- £ (e), 06) 

and the bottom loss per bounce can be defined as 

« w (0 = -201og(j*(e „))= 20e(e„). (37) 

Combining the equations above, the reduction in the reflection coefficient amplitude can be 
determined by 

' 20 iV(0„) 20MAf(e„)' ( ) 

Bottom attenuation can be introduced by defining an imaginary component of the 
bottom sound speed, 



c b-* c b~ & > (39) 

such that the wavenumber in the lower medium becomes 

k b =-->k b +ia e , (40) 

c b 

where 



a =- 



©5 



(41) 



Since the solution in the lower medium goes like 



e ikhR -» e k * R e a ‘ M . 



(42) 



then a e is found to produce the exponential decay in range. This then causes transmission 
loss according to 



A TL atten = 8.686a, A R = a Aft , (43) 

where a is the bottom attenuation in units dB/m, and can be written as 

a = 8.686^. (44) 

C b 

In order to estimate the bottom attenuation (dB/m), a simple matching algorithm 
may be employed which searches over values of 6 (the complex bottom sound speed) until 
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the reduction in the reflection coefficient amplitude agrees with the observed value of e(0 n ) 
defined in Eq. 37. Once values for 8 are obtained (for various frequencies and mode angles 
9 n ), then Eq. 43 may be employed to produce the final estimate for the bottom attenuation 



a. 
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IV. DATA ANALYSIS AND RESULTS 



A. INVERSION OF BOTTOM SOUND SPEED 

As described in Section III. A., the inversion for bottom sound speed is based on the 
evaluation of the vertical coherence of the ambient noise field. For this analysis, the 
ambient noise was extracted from segments of the SUS data sets between actual SUS 
signals (and allowing for reverberation to fade). The spectrum of the ambient noise 
computed showed little variation from low frequency up to 600 Hz indicting no strong 
unexpected noise sources were present. Also, the imaginary component of the coherence 
between the two hydrophones was found to be small as required. This indicates that the 
ambient noise field at the trial site for frequencies below 600 Hz has a symmetrical 
distribution in the vertical direction and the data is suitable for inversion by the method 
described previously. 

It should be noted that the inversion method described in Section III. A. was 
originally defined for an isospeed water column overlying a homogenous bottom. As is 
obvious from the local sound speed data (refer to Figure 2), there is significant variation in 
sound speed over depth. However, if the sound speed change is not too large over the 
vertical separation of the phones, the inversion method should still be applicable by simply 
adjusting the local angles to the appropriate values at the bottom interface using Snell’s law. 

Since the data from the upper two phones of the VLA was analyzed, a local sound 
speed of c w w 1485 m/s was used in the calculation of the ambiguity function (refer to Eq. 
11 ). 

The fit between the theoretical coherence obtained from the model and the real 
component of the observed coherence is shown in Figure 7. By computing the ambiguity 
surface for the optimization of two parameters (a 2 and 0 f ) an optimal value was found at 

0 C » 1 8° in the water column. Figure 8 shows the ambiguity surface for the optimization 
with two parameters, a 2 and 0 C whose minimum corresponds to a 2 « 0. 1 8 and mean critical 
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angle 0 f *18° ± 1° for 4 different SUS explosive data. A grazing critical angle of 
approximately 19°±1° near the bottom was determined from Snell’s law. The sound speed 
close to the bottom was 1476 m/s, and therefore the bottom sound speed c 2 «1563±10 m/s. 
This value was compared to values obtained by a sophisticated genetic algorithm technique 
employed at the University of Rhode Island (Miller and Potty, 1998). They found values in 
the upper few meters of the sediment in agreement with the result reported here. 

B. INVERSION OF BOTTOM DENSITY 

To invert for bottom density, the method of ideal waveguide effective depth was 
introduced. However, this method assumed a constant sound speed in the water column. In 
fact, this was fundamental in the definition of the ideal waveguide eigenfunctions. For that 
reason, this method may not be expected to produce very good results. We will attempt to 
invert for density anyway by assuming the water column is isospeed with nominal sound 
speed c w « 1485 m/s. Consistent with this value, we shall assume the critical angle 0 C « 1 8°. 
In order to obtain separated mode travel times, spectrograms were computed for various 
SUS explosive data. The theoretical mode travel times, for the ideal waveguide were then 
computed as a function of range R between source (SUS) and receiver (VLA), water depth 
H, water sound speed c w and critical angle 0 C . A comparison between theoretical and 
measured travel times for all 1 1 SUS data sets produced an average value for the density 
ratio of 1 .6. However, the changes in the theoretical curves for various density ratios was 
not large whereas the mode arrival structures were fairly ambiguous. Thus the uncertainty 
is very large. It is also possible that the theoretical mode travel times have significant errors 
due to the mismatch in sound speed profile. However, the estimated value of p b « 1.6 g/cm 3 
from the inversion method seems quite realistic. The results supplied by Miller and Potty 
suggest a more reasonable estimate for the bottom density is about 1.8 g/cm 3 . Figures 9a 
and 9b show the mode travel times for a phone from a SUS explosive data and those for an 
ideal waveguide. 
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C. INVERSION OF BOTTOM ATTENUATION 

A method was developed in Section III.C. to invert for bottom attenuation using 
measures of modal transmission loss as a function of range. In order for this method to 
work successfully, the decrease in amplitude due to bottom attenuation must be observable 
over the transmission loss due to cylindrical spreading. For example, as the range between 
the SUS drop location and the VLA increases from 15 km to 30 km, the change in TL due 
to cylindrical spreading is approximately 3 dB. This equates to about 0.2 dB/km. In order 
for bottom attenuation to be distinguishable, Eq. (38) suggests that the quantity 20 N s 
should be of this same order of magnitude or larger. 

A simple analysis assuming bottom attenuation increase linearly with frequency 
shows that s increases with increasing frequency. Specific values can be obtained using 
Eqs. (34) - (44) with realistic values of bottom attenuation a. However, as frequency 
increases, the angle of propagation for a specific mode decreases thereby decreasing N, the 
number of bottom interactions per km. 

If one assumes a reasonable value of a ~ 0.05 dB/km/Hz, a 100 Hz signal suffers a 
0.5 dB/km loss due to attenuation for angles of propagation of about 12°. At this frequency 
in this waveguide, that corresponds roughly to Mode 2. Unfortunately, this mode is not 
clearly distinguishable at this frequency from the available SUS data. The lowest modes 
can be isolated at lower frequencies, but then the attenuation is negligible. At higher 
frequencies, the propagation angles for these low modes is near grazing, and again the 
influence of attenuation is small. 

The result of this analysis, unfortunately, shows that this data does not avail itself to 
inversion for bottom attenuation using the method developed previously. Other processing 
methods may be able to extract the necessary information in order to utilize this method, but 
will not be explored further here. 
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Figure 7. Coherence of the Measured Ambient Noise between Two Hydrophones; Real 
Component (Dashed Curve) and Imaginary Component (Dotted Curve). A Good Fit 
Between the Theoretical Coherence Obtained from the Model (Solid Line) and the Real of 

the Measured Coherence is Observed. 
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Figure 8. Ambiguity Surface for the Optimization with Two Parameters, a 2 and 0 C . The 
Minimum is Marked by X which Corresponds to 0.18 and 0 C « 18° for 4 Different SUS 

Explosive Data. 
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Figure 9a, Mode Travel Times for a Phone from a SUS Explosive Data 




Figure 9b. Mode Travel Times for a Phone from an Ideal Waveguide. 
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V. CONCLUSION 



This thesis incorporated several methods which may be applied directly to measured 
data to provide a relatively simple means of obtaining reasonable estimates of ocean bottom 
parameters from minimal information. The analysis of the data from the Mid-Atlantic 
Bight experiment recorded on the NW VLA produced the following results: 

1 . Ambient noise was used to measure both the vertical directionality and spatial 
coherence of the noise field. By examining the noise data, we have been able to 
estimate 0 C and hence the sound speed c b in the sediment. The bottom sound 
speed deduced was approximately 1563110 m/s. 

2. The method of effective depth was employed to estimate the bottom density p b . 
Although this method was derived for an isospeed water column environment, it 
is expected to produce reasonable estimates for most shallow water 
environments. By comparing mode travel times extracted from SUS data 
spectrograms and comparing these with values for an ideal waveguide with an 
effective depth, the bottom density was estimated to be p b « 1.6 g/cm 3 . The 
uncertainty in this value is very large, however, partly due to the difficulty in 
distinguishing mode arrivals clearly in the spectrograms. 

3. In order to distinguish bottom attenuation from simple cylindrical spreading, it 
was determined that values for mode transmission loss were required for high 
modes at high frequencies. Unfortunately, the SUS spectrogram data could not 
adequately provide this information, and no estimate of bottom attenuation was 
produced. It may be possible to use the method outlined if a modal 
decomposition of the data on the VLA could be performed. However, due to the 
relatively short aperture of the array (compared to the waveguide depth), this 
seems unlikely. Furthermore, the source level of different SUS detonations may 
not be repeatable enough to use this data for accurate estimations of bottom 
attenuation. 
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Comparing the inversion results for bottom sound speed and density with values 
obtained from genetic algorithm techniques provided by the University of Rhode Island, 
good agreement was found particularly for the sound speed. However, the analysis of the 
SUS data proved more difficult than originally anticipated. Much of this has to do with the 
uncertainty in absolute drop location and repeatable source signals. Another part of the 
problem was the reliance on low mode information which is difficult to characterize for 
SUS data measured on only eight hydrophones with total aperture less than half the width 
of the waveguide. It is likely that this analysis would be more successful with different 
more well-defined sources and a full water column array. 
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